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A new model that generalizes the study of quantum Brownian motion (BM) is constructed. We 
consider disordered environment that may be either static (quenched), noisy or dynamical. The 
Zwanzig-Caldeira-Leggett BM-model constitutes formally a special case where the disorder auto- 
correlation length is taken to be infinite. Alternatively, localization problem is obtained if the noise 
auto-correlation time is taken to be infinite. Also the general case of weak nonlinear coupling to 
thermal, possibly chaotic bath is handled by the same formalism. A general, Feynman- Vernon 
type path-integral expression for the propagator is introduced. Wigner transformed version of this 
expression is utilized in order to facilitate comparison with the classical limit. It is demonstrated 
that non-stochastic genuine quantal manifestations are associated with the new model. It is clarified 
that such effects are absent in the standard BM model, either the disorder or the chaotic nature 
of the bath are essential. Quantal correction to the classical diffusive behavior is found even in the 
limit of high temperatures. The suppression of interference due to dephasing is discussed, leading to 
ff) ■ the observation that due to the disorder the decay of coherence is exponential in time, and no longer 

depends on geometrical considerations. Fascinating non-Markovian effects due to time-correlated 
(colored) noise are explored. For this, a new strategy is developed in order to handle the integration 
over paths. This strategy is extended in order to demonstrate how localization comes out from the 
| path integral expression. 

I. INTRODUCTION 

The dynamics of a particle that interacts with its environment constitutes a basic problem in physics. Classically, 
upon elimination of the environmental degrees of freedom, the reduced dynamics is most simply described in terms of 
Langevin equation. Solution of this equation, by utilizing Fokker Planck equation, is well known JO]. In the absence of 
external potential, it yields spreading and diffusion. The latter effect is due to the interplay of noise and dissipation. 
However, diffusion may arise also from the interaction with disordered environment. This kind of non-dissipative 
"random-walk" diffusion is encountered, for example, in Solid State Physics, while analyzing electrical conductivity. 
It is well known ]ll| that this latter type of diffusion may be suppressed quantum mechanically due to localization 
■ effect. Still, diffusive- like behavior is recovered if noise and dissipation are taken into account. 

The unified modeling of the environment in terms of noise, dissipation and disorder is the first stage of the present 
. study. One may take the notion of particle literally, and identify the environment as cither external or internal bath 
that consists of infinitely many degrees of freedom. The bath may be a large collection of other particles or field 
modes (photons, phonons). Else, the internal degrees of freedom of the particle itself are considered to be the bath. 
The latter point of view has been suggested by Gross || in order to analyze inelastic scattering of heavy ions. 

A totally different point of view, promoted by Caldeira and Leggett (CL) |t| , consider the notion of particle as a 
token for some macroscopic degree of freedom. A linear interaction with a speculated bath that consists of infinitely 
many uncoupled harmonic oscillators, is assumed (Zwanzig ||). The known classical limit, namely Langevin-type 
equation, serves as a guide for the construction of the proper Hamiltonian. (Phcnomcnological rather than microscopic 
considerations are used, hence the usage of the term 'speculated bath'). The power inherent in this approach is the 
capability to introduce an explicit path integral expression for the propagator, using the Feynman Vernon (FV) 
|^| formalism. This FV-CL propagator constitutes a quantized description of Brownian motion (BM). The term 
"BM-model" will be associated from now on with this propagator. 

The first question which should be asked concerning the applicability of CL approach is obviously whether either 
the coupling, the bath, or both are "too simple" in order to account for the variety of physical phenomena that are 
associated with generalized BM. The term "generalized BM" is used in order to describe any dynamical behavior that 
corresponds in the classical limit to Langevin-like equation. In the simplest case Langevin equation is mx + rjx = T 
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where m and r\ are the mass and the friction coefficient, respectively. One should specify the stochastic force T . This 
force may arise from interaction with some fluctuating homogeneous field (T — !F{ij), which is the usual formulation. 
However, more generally, this force may arise from the interaction with disordered potential (T = — VW(x, t)). In the 
latter case the spatial auto-correlations of the force are significant. To avoid misunderstanding it should be emphasized 
that there are other aspects in which BM can be generalized (for review see M). 

In the present paper we construct a unified model for the study of Diffusion, Localization and Dissipation (DLD). 
This model describes generalized BM in the sense specified above. The disordered environment may be either static 
(quenched), noisy or dynamical. The model is treated within the framework of the FV formalism. The resultant path 
integral expression for the propagator contains a functional Sf with kernel a(r— r') that corresponds to friction, and a 
functional S n with kernel 4>(t— t') that corresponds to the noise. Both functionals depend also on a suitably defined 
auto-correlation function w(r—r') that characterizes the disorder. BM model constitutes (formally) a special case 
where the disorder auto-correlation length is taken to be infinite. Optional derivations of the resultant path integral 
expression are presented for the particular cases of either classical or quantal systems with non-dynamical disordered 
potential. Obviously, in the latter case the friction functional Sf is not generated by the derivation. Localization 
problem is obtained if the noise auto-correlation time is taken to be infinite (4>(t— t') = const). 

An explicit computation of both the classical and the quantal propagators )C(R, P\Rq, Pq) will be carried out. This 
propagator generates the time evolution of p(R, P), which is either the Wigner function or the corresponding classical 
phase-space distribution. Spreading and diffusion profiles are found for either noisy or ohmic environment. Quantal 
corrections to the classical result are discussed. A new strategy is developed in order to handle the integration over 
paths. This strategy is utilized in order to study the anomalous "diffusion" profiles due to colored noise. Later it is 
extended in order to demonstrate how localization comes out from the path integral expression. 

Again, one may ask, whether the DLD model is the "ultimate" model for the description of BM in the most 
generalized way (as far as generic effects are concerned). In case of 2-D generalized BM one should consider also the 
effect of "geometric magnetism" pij , which is not covered by the 1-D DLD model. Here we limit the discussion to 1-D 
BM. In order to answer this question one should consider a general nonlinear coupling to a thermal, possibly chaotic 
bath. In the limit of weak coupling one may demonstrate (see Sec.VI(B)) that indeed the bath can be replaced by 
an equivalent "effective bath" that consists of harmonic oscillators, yielding the DLD model. A further reduction to 
BM model is achieved if the coupling is linear. The derivation also demonstrates why the so called "ohmic" bath is 
generic. However, we cannot prove that the path-integral expression that corresponds to the DLD model is the most 
generalized description of BM (in the sense of this paper). Gefen and Thouless ||, Wilkinson 10 and Shimshoni and 
Gefen || have emphasized the significance of Landau-Zener transitions as a mechanism for dissipation. The weak 
coupling approximation misses this effect. Still, there is a possibility that some future, more sophisticated derivation, 
will demonstrate that an equivalent "oscillators bath" can be defined also in the case of strong coupling. The existence 
of such derivation is most significant, since it implies that no "new effects" (such as "geometric magnetism" in case 
of 2-D generalized BM) can be found in the context of 1-D generalized BM. Referring again to the Landau-Zener 
mechanism, Wilkinson has demonstrated that anomalous friction, which is not proportional to velocity, may arise 
fLOf . The BM model cannot generate such anomalous effect. However, we shall demonstrate that the non- ohmic DLD 
model can be used in order to generate this effect. 

We turn to review previous works that are related to the present study. Traditional approach to the study of 
dissipation is based on the Master Equation formalism B , which is the quantal analog of the classical Fokker Planck 
treatment jD. Systematic derivations are usually based, in some stage, on the Markovian approximation. An alterna- 
tive route is to apply the FV formalism || . This has been done by Mohring and Smilansky | (| , following Gross || , in 
order to study deep inelastic collisions of heavy ions. Later, the FV formalism has been applied by CL [Q and followers 
Q in particular to the study of macroscopic quantum tunneling. Hakim and Ambegaokar Q| has applied FV model 
in order to compute the spreading and the diffusion of quantum Brownian particle. Cohen and Fishman pl[ | have 
computed the full Wigner propagator in case of general quadratic Hamiltonian, possibly time dependent. The latter 
study has demonstrated more clearly the significance of noise time- autocorrelations. Non-trivial auto-correlations 
may lead to non-Markovian effects due to the non-local (in time) nature of the noise functional Sn- The simplest 
example is the suppression of diffusion at low temperatures |8) due to negative power-law correlations of the noise 
fHf . A less trivial manifestation of non-Markovian effect has been found by Cohen |23| while analyzing diffusion due 
to the destruction of localization by colored noise. 

Non Markovian features are usually speculated to be less relevant if disorder is taken into account. One may 
expect that for generalized BM, disorder will lead, at all temperatures, to normal diffusion. We shall find later in this 
paper that non-Markovian effects are quite effective also in case of the DLD model, and lead to anomalous spreading 
profiles. However, it will be clear that these effects, though counter intuitive at first sight, are of classical nature. 
Non Markovian effects can also arise from the non- locality of the friction functional Sf, this is the case for non-ohmic 
bath. The retarded response of the bath has then long memory for the particle's dynamics, see for example the review 
papers by Grabert et al. and by Hanggi et al. 0]. In particular, for unbounded motion, in the absence of disorder 
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(BM model), it will be demonstrated that one encounters infinities in computations of the friction and of the effective 
mass. These infinities are avoided if disorder is taken into account (Sec.IV(B)). 

The suppression of quantal interference due to dephasing process is an important issue for both semiclassical |^] 
and mesoscopic |Tg| l physics. The dephasing of interference in metals due to electromagnetic fluctuations has been 
discussed by Al'tshuler, Aronov and Khmelnitskii ]l7|] . General considerations has been presented by Stern, Aharonov 
and Imry [Q. The strong dimensionality dependence of the dephasing process has been emphasized. In this paper, 
the DLD model is used in order to study the suppression of quantal interference due to the local interaction with either 
noisy or dynamical disordered environment. The dephasing is determined by the noise functional Sn[x"(t), x'(t)]. 
Due to the local nature of the dephasing process, the decay of coherence is exponential in time, and no longer depends 
on geometrical considerations. 

In case of quenched disordered environment the DLD model reduces to a localization problem that its solution is 
well known [ 12]] . In particular, for 5 correlated potential the result for the localization length has been pointed out by 
Thouless ]TJ. As far as we know, the real-time Feynman path integral formalism has not been utilized so far in order 
to re-derive this result, though functional integration is frequently employed in closely related computations. The case 
of noisy disordered environment, where the potential is 6 correlated both in time and in space has been considered 
by Jayannavar and Kumar p5| . In the latter reference the spatial spreading has been computed, and a classical- like 
result has been obtained. Quantal corrections to the dispersion profile, has not been discussed in the latter reference. 
No solution exists for a model that "interpolates" the crossover from quenched noise, via colored noise, to white noise 
disordered potential. Marianer and Deutsch |f6| have considered the problem of white noisy disordered potential with 
added dissipation. Using BM model, they have demonstrated that a classical-like results is obtained for the spatial 
spreading. Again, neither the spreading profile nor quantal corrections have been considered. 

The outline of this paper is as follows: In Sec. II the unified model for the study of Diffusion, Localization and 
Dissipation (DLD) is constructed. Derivation of the reduced classical dynamics is presented in Sec.III-IV. It is shown 
that a well defined Langevin equation is obtained for both subohmic and superohmic bath, as well as for ohmic 
bath. This is contrasted with BM model, where only the ohmic case is well defined. In Sec. V- VI, four derivations 
of the FV path integral expression for the propagator are presented: (a) A classical derivation that is based on 
Langevin equation; (b) A quantum mechanical derivation for non-dynamical noisy or quenched disordered potential; 
(c) A quantum mechanical derivation for the full DLD model; (d) A quantal derivation for the general case of weak 
nonlinear coupling to a thermal, possibly chaotic bath. App.A clarifies the relation of BM-model and DLD-model to 
the case of interaction with external bath that consists of extended field modes. In Sec.VI(C) the explicit expression 
for the influence functional allows concrete predictions concerning the loss of interference due to dephasing. App.B 
introduce a gedanken experiment that clarifies the manifestation of interference and dephasing in the presence of 
dynamical disordered environment. In Sec. VII- VIII three strategies for the computation of the quantal propagator 
are discussed. Spreading and diffusion profiles are found for either noisy or ohmic environment. The DLD model 
is compared with the BM model. In the latter case the result is classical-like, while in case of the DLD model, a 
singular "quantal correction" should be included. The significance of this "correction" is further clarified in App.B. 
The strategy of Sec. VIII, which enables computation of diffusion profiles in the presence of colored noise, is extended 
in Sec. IX. We use this strategy in order to demonstrate how, in the case of quenched noise, localization comes out 
from the general FV path integral expression. Summary and Conclusions are presented in Sec.X. 



II. CONSTRUCTION OF THE DLD MODEL 



A. Langevin Equation 

As a starting point for later generalization we consider the classical Langevin equation. This equation describes 
the time evolution of Brownian particle under the the influence of so called "ohmic" friction and stochastic force. 

mx + rjx + V (x) — J- [Langevin Equation] . (2-1) 

In the standard Langevin equation the stochastic force represents stationary "noise" which is zero upon ensemble 
average, and whose autocorrelation function is 

(T{t)T{t')) = 4>{t-t') [for standard Langevin]. (2.2) 

Usually white noise, which is 5 correlated in time, is considered. The standard Langevin equation can be generalized 
by assuming that the stochastic force is due to some noisy potential, namely T{x,t) = —U'(x,t) where / denotes 
spatial derivative. One may assume that U is zero upon ensemble averaging, and satisfies 

(U(x,t)U(x' ,t')) — (f>(x— x',t— t 1 ) [generalized Langevin] (2.3) 
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The autocorrelation function of the stochastic force T at some specified point x is 4>(t—t') = — <fi"(0, t—t') In practice 
cj>(x—x',t—t') will be assumed to have the factorized form 

4>{x-x\ t-if) = (f>(t-t') ■ w{x-x') (2.4) 

where, without loss of generality, w(r) is normalized so that w"(0) = —1, the normalization constant being absorbed 
into 4>{t). 

For distribution of particles, the standard Langevin equation predicts rigid motion with no diffusion. This feature 
is eliminated if an average over realiza tion s of T is performed. The standard Langevin equation can be viewed as a 
special case of the generalized version (2.3). For this one should take hi{x,t) — —x ■ T(t). Alternatively, 

w(x — x') = — — (x — x 1 ) 2 [for standard Langevin]. (2-5) 



Another choice for spatial autocorrelations that corresponds to disordered environment is 

w(x — x 1 ) = I 2 exp ( — — ( — -,, — ) ] [for disordered environment]. (2-6) 



The parameter I constitutes a measure for the microscopic scale of the disorder. For distribution of particles, the 
resultant motion will be diffusive-like rather than rigid, even without averaging over realizations. It is crucial to do 
this important observation if one wishes to introduce a quantized version of Langevin equation. 



B. Hamiltonian formulation 



Disregarding the friction term, the Langevin equation may be derived from the Hamiltonian 
V 2 

H=^- + V(x)+U(x,t) (2.7) 

where average over realizations of U is implicit. If hi is time independent, we shall use the notion quenched disordered 
environment. If hi is uncorrelated also in time, we shall use the notion noisy environment. 

If rj is non-zero, energy may dissipate, and we should consider "dynamical" environment rather than "noisy" or 
"quenched" one. In order to model the environment we generalize the FV-CL approach. Namely, to facilitate later 
mathematical treatment, the interaction potential is assumed to be linear in the bath coordinates: 

U(x,t) -> Hj nt = -^c a Q a u a (x) (2.8) 

a 

where Q a denotes the dynamical coordinate of the a scatterer or bath mode. u a (x) is the normalized interaction 
potential, and c a are coupling constants. In case of dynamical environment, the Q a are assumed to be oscillators' 
coordinates. The bath Hamiltonian is: 

= + \ m "lQD ■ ( 2 - 9 ) 

a a 

The total Hamiltonian of the system plus the heat bath is 
V 2 

n = ^- + v(x) + n int + n bath . 2.10 

In order to further specify the DLD model, we should characterize the spectral distribution of the bath-oscillators, as 
well as their interaction with the particle j2q] . 

In what follows we consider the case of localized scatterers. The scatterers are assumed to be uniformly distributed 
in space with u a (x) = u(x—x a ). The joint distribution of the bath oscillators with respect to uj a and x a is assumed 
to be factorized. Consequently we can write 

2 

~ -^—8(U) - u a ) S(X - x a ) = J{u) . (2.11) 
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The bath is characterized by the spectral function J(oj). If we consider a partition of the bath oscillators into subsets 
of oscillators whose positions x a are the same, then locally, within each subset, the uj a will have the same distribution. 
The interaction is characterized by well defined spatial autocorrelation function w(x-x'), namely 

/oo ^ 
u(r - x')u{x')dx' = density z2 u <x{R + r)u a (R) , (2.12) 

where "density" refers to the uniform spatial distribution S(x—x a ). The scattering potential u(x) will be normalized 
so that w"(0) = — 1, the normalization constant being absorbed in the coefficients c a . Disregarding for a moment 
the dynamical nature of the bath, thus considering again the case of either "noisy" or "quenched" environment, one 



obtains (U(x,t)U(x r ,t')) = ^ a c a {Q a (t)Q a (t'))u a {x)u a {x') which upon recalling previous definitions, leads to (2.4) 
with cf>(t—t') = density x J2a c a(Qa(t)Qu(t')) ■ Thus, if the dynamical nature of the bath is ignored, the problem 
reduces to solving Langevin equation with the appropriate w{x— x') and <fi(t—t'). 

The spectral functi on J (u>) and the spatial autocorrelation function w(r) constitute a complete specification of the 
model Hamiltonian ( |2.10| ). This Hamiltonian will be used in order to study diffusion localization and dissipation 
(DLD) within the framework of a unified formalism. 

The BM model, where u a (x) = x up to normalization, constitutes formally a special case of the DLD model. It 
will be possible to apply the "unified" formalism that will be developed for the DLD model, al so f or the analysis of 



the BM model, merely by substitution of the appropriate w(r). Namely, one may use w(r) of (2.5), or alternatively 
one may take the limit I — > oo in (^^). However, there are some minor, but significant differences that will be noted 
in due time. This is because our derivations rely on the local nature of w(r). Therefore we consider the notion "BM 
model" distinctively from the notion "DLD model" . Also the results will be both qualitatively and quantitatively 
different, in-spite of being derived, formally, from look-alike formulas. 

Fig.|l| illustrates the BM Hamiltonian and the DLD Hamiltonian. Looking at the figure it is apparent that the 
natural application of the BM model is for the description of the dynamics of a composite particle, with many 
internal degrees of freedom. This interpretation is due to Gross Ref. [|| , who suggested to apply such models for the 
analysis of heavy-ion collisions. The bath according to this interpretation is "internal", carried by the particle, rather 
than "external" . 



The interaction Hamiltonian (2.?) may also describe the interaction with some "external" bath that consists of 
extended field mo des . For example, the electron-photon interaction and the electron-phonon interaction may be cast 
into the form of ( |2.8| ), with u a (x) = cos(u) a /c ■ x), sin(u> a /c ■ x) , where c is either the speed of light or the speed 
of sound. In the long wavelength limit this interaction resembles the BM-model rather than the DLD-model. With 
minor modifications, also these models may be treated within the framework of our "unified" treatment (Appendix 
A)-_ 

Finally, we re-emphasize that in general, the DLD model as well as the BM model may be used on a phenomenological 
basis for the description of dissipation in mesoscopic quantum devices. This point of view will be discussed further 
later on (Sec.VI(B) in particular). 



III. DERIVATION OF LANGEVIN EQUATION 

A. The Reduced Equation of motion 

The quantal state of the particle may be represented by Wigner function p{R, P). The time evolution of Wigner 
function corresponds to that of classical distribution in phase-space. In this section we consider a 'classical treatment 
of the dynamics'. The latter term implies that the system is considered to be classical, (h sys t e m 0), while the bath 
gets full quantum mechanical treatment. The limit htath ~ ► is not taken. The equations of motion of classical points 
that form a distribution in phase-space are x = p and 

P = -V {x) + T ha th , Fbath = ^C a Q a (t)u' a (x) . (3.1) 

a 

The variables Q a {t) satisfy the equation 
m a Q a {t) + m a ujlQ a (t) = c a u a (x) 
which can be solved explicitly, namely, 

Qa{t) — Qa{0) cos(uj a t) H sm(ui a t) + / dt' — - — sin(w Q (i — t'))u a (x{t')) . 

m a uj a J Q m a uj a 
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Substitution of the latter expression into ( p.l| ) yields 

-Pbath — retarded ~t~ P noise 

with 

^retarded = / dt' 2a{t~t') w' (x (t) - X (t' )) 



The response kernel a(r) is defined for positive times (0 < t) as follows: 

dto J(lo) 



a(r) 



du) d 
— J (ui) sm(u!T) = — — 
7T or 



■ cos(o;r) 



(3.2) 



(3.3) 



(3.4) 



In order to make further progress, a specification of the initial state of the system plus the bath is needed. We shall 
assume that initially (at time t = 0) the system is prepared in some arbitrary quantal state while the bath oscillators 
are in thermal canonical equilibrium with some reciprocal temperature f3. The Wigner function representation of the 
probability density matrix is then (see App. A of Ref. fl22j): 



Pi 



=0 (R, P;Q a ,P a ) = Pt=o(R, P) • \{ p eq {Q a , P a ) 



where 



Peq (Qa j Pa 



■ exp 







/tanh(i^/3w Q )\ / P 2 



i coth(A?i/3a;, 
Using (|3.6[) one obtains the expectation values 



-m a ujlQl 



\2m a 



(3.5) 



(3.6) 



P (0) 2 1 11 

= (o m «^Q a (0) 2 ) = -hu; a coth(-hf3La a ) . 
2m a 2 4 2 

Hence it is easily found that {J~ no i se {t)) = while 

{T notse {t)T notsel {t')) = -J2cl(Q a {0) 2 )cos(uj a {t-t'))u' a {x)u' a (x') . 

Oi 

Thus, one identifies that the noise term is characterized by the autocorrelation function ( |2.4| ) with 

l>OG 

<j>{t-t')= I —J(uj)hcoth(-h/3u) cos[w(i-i')] 

./n 7T 2 



(3.7) 



The Langevin equation (3.1) together with ( 3.2 )-( 3.4) and ( |3.7| ) constitutes an exact and complete description of the 
reduced dynamical behavior of the system, as long as the system is considered to be classical in nature. 



B. Langevin Equation in the Absence of Disorder 



For the BM model u a (x) = x. The derivation of L angevin equation (|2 
The latter may be interpreted as a special case of (2.3), provided equation (12 
expression for the retarded force (3.3), w'(x{t)— x(t')) is replaced by x(t'), rather 
significance of this minor difference is discussed below. In the equation for the 
the integration over t' to infinity, provided a(r) is replaced by <S(t) = a(r)9(r) 
turn, this kernel may be written as a sum of three terms, namely <5(r) = ao(r) 
the asymmetric continuation of a(r) to the domain r < 0, while ao(r) + a m (r) 
latter is split into ao(r) which is a delta function, and a m (T) which satisfies J 
retarded force, for the BM model, may be written as the sum: 



leads to T n oise that satisfies fl2.2p . 

t) for w(r) is used. However, in the 
than by — (x(t)— x(t')). The physical 
retarded force (3.3) one may extend 
, where O(r) is the step function. In 
+ a(r) + a m {r) . The kernel a(r) is 

is the symmetric continuation. The 
a m (r)dT = 0. Consequently, the 



■P retarded — -P switching "f" P. {^potential ~t~ *P friction ~t~ P, Amass \^*^J 

where F SW itching ~ -k -6(t)-x(t), and Potential ~ +Ak-x(t), and J 7 friction ~ -r) e ff±(t), a,ndJ 7 Amass -Am-i(t), 
with fco = ?7 e // = lim^^o an d Afc = ^ '■^^-duj 1 and Am = — a{r) r 2 dr. It has been assumed that a(r) 
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has short range duration r c , much shorter than the physically relevant time scales of the dynamics. This assumption 
is not true in general, and the consequences will be discussed later in the next section. 

The switching impulse act on the particle if it starts its trajectory at a point x ^ 0. This term originates due to 
the fact that the initial preparation is such that the bath is in thermal equilibrium provided x — 0. The ^potential 
force may be avoided if we care to include in the BM Hamiltonian the proper "counter term". Namely: 

n bath + n mt -> E(^ + H^-^ 2 ) (3J) 

This expression is manifestly invariant under space translations. Sanchez-Canizares and Sols have further considered 
this issue p5[ . 



C. Langevin Equation for Disordered Environment 



We turn back to the DLD model with J- re tarded as found in equation (3.3). Here, performing similar treatment, 



neither ^switching nor 3~ Apotential are encountered. Formally, this is due to the fact that the difference —(x(t)—x(t')) 



appears in (3.3), r ath er than x(t') by itself. Physically, this is due to the inherent homogeneity of the interaction with 



the environment ( |2.8| ). As for J 7 friction and T Amass, here one obtains 

/oo 
orfr) w'(vt) dr 
-OO 

a(r) w"(vt) t 2 dr (3.10) 



If a(r) has short range duration r c , and if furthermore the velocity of the particle is not too hig h (v ■ t c « <), then, 
using w'(0) = and w"(0) = —1, the above expressions reduce to those that follow equation ( |3.8| ). If a(r) is not 
short range, to be discussed below, then the retarded force in the BM model is characterized by long range memory 
for the dynamics. The approximation that has been used in the preceding subsection is no longer valid, and infinities 
are encounter ed in the co mpu tations of the constants there. In general, these annoying features are not shared by the 
DLD model ( fl3TC| ) with (|J)), since a finite cutoff T£ = t/v exists, where the microscopic length scale t characterizes 



the interaction range with the scatterers. 



IV. OHMIC AND NON-OHMIC BATHS 



In o rder to make further progress, a specification of the spectral function J(w) is required. This function, defined in 



( 2.11 ), characterizes the distribution of the bath-oscillators with respect to their frequencies. Following CL we assume 



it to be of the form 

J(ui) — t]lu s G(uj/lo c ). for "subohmic" bath < s < 1 

for "ohmic" bath s = 1 

for "superohmic" bath 1 < s (4-1) 

Here the exponent s characterize the singular behavior of J{lo) in the vicinity of lj = 0, while G denotes a smooth 
cutoff function. The latter satisfies G(0) = 1. Moreover, it is assumed that G(|w|) is analytic for real frequencies. 
For example, G may be chosen to be either Lorentzian or Gaussian. The latter possibility will be adopted from now 
on. The commonly used Exponential cutoff exp(—Lo/oj c ) will not be used since it does not satisfies the mentioned 
requirements. Exponential cutoff results in singular behavior that corresponds to the exponents s, (s + 1), (s + 2), ... 
rather than "pure" singular behavior that corresponds to s alone. The asymptotic behavior of both a(r) and 4>(t) 
is dictated by the singularity of J{oj) at lu = 0. The physical significance of s is further discussed in App.A and in 
Sec.VI(B). 

A. Expressions for the Kernels 

In order to display explicit expressions for both a(r) and </>(t), it is useful to define the following function: 

/•OO 7 

G s (t) = / — u 3 - 1 e~^ )2 cos(wr) (4.2) 
Jo n 
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1 2 3 4 5 6 

time [normalized units] 

FIG. 2. Plots of G s (t) for s = 0.5, s = 1, s = 2 and s = 4. The plots are easily distinguished by referring to their description in the 
text. 



For s = 1 this function is a normalized Gaussian. For odd s one obtains 

d 2n 



G 1+2n (r) = (-1)' 



[n is integer] 



For general s it starts at r = positive 

g s (o) = i- 2 s / 2 r(~) ^ 



(4.3) 



(4.4) 



and cross the zero value [(s + l)/2] times (now we consider s ^ odd number, possibly non- integer) . The short range 
oscillatory behavior dies out on a time scale of the order r c = a;" 1 . For larger times a power-law decay is found: 

G„(t) = (- cos( Js)r(s) ) • — [tail for s ^ odd] 



(4.5) 

The total algebraic "area" under G s (t) is infinite for < s < 1, finite for s = 1, and zero for 1 < s. Representative 
plots of G s (t) are displayed in FigJ|. 

For the spectral function J(uj) as in (4.1), with Gaussian cutoff, the kernel a(r) is 



a(r) = ~rj -^G s (t) 



(4.6) 



Also 4>(t) can be expressed in terms of G s (t) in both cases of "high" and "low" temperatures. In the first case it is 
assumed that Hf3 (that has dimensions of time) is much shorter than any dynamical time scale. One may use then 
the approximation 



(t) = 2^ G s (t) ["high" temperatures]. 



Else, if Hf3 is long, one obtains 
0(r) 



|r (s+ i) C ( s+ i) • ^ (i 



s+l 



TirjG s +i(T) ["low" temperatures], 



(4.7) 



(4.8) 



where £(s) is Riemann Zeta function (sum over l/n s ). The notion "low temperatures" means here r < h/3. The 
temperature-dependent constant term results from the integration ^ J(u})h[coth.(^hf3uj) — 1]. From now on we 
shall use the notions of "high" and "low" temperatures in the sense of the above approximations. 
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B. Expressions for Friction and Effective Mass 



We turn back to the computation of the terms in (3.8). For the Apotential term the result is always finite 



1 s 

Ak = -2 s / 2 T(-) nuj s c [switching impulse in BM model] 



(4.9) 



For the friction one obtains a finite, non-zero result, (r) e ff — rj) only in the case of an ohmic bath (s = 1). For 
subohmic bath (0 < s < 1) the friction is infinite, while for superohmic bath (1 < s) the friction is zero. In case of 
the DLD model finite results are obtained for all cases. Using ( fUOD one finds out 



J~ friction {p 

2 



ryv 



dr 



du> 



d r 



— u s e 2U <= ; shhwr) • — 
7r ar 



T ( C 



2 \ t. I 



rn+-i 2 s/2 ± 



-r(i+-) 2 s / 2 

7T 2 



nv 



(4.10) 



Above, the notation Tg = £/v has been used. In the last step the limit tu c — > oo has been taken since it leads to a 
finite, non-zero result. Note that the cutoff frequency uj c is not relevant physically as long as v -C £ ■ u) c . The above 
computation clarifies the significance of the various time scales. Formally, the BM model constitutes a special case 
with T£ = oo. For s = 1, the so called "ohmic" case, one obtains !Ff r iction — —nv. This result holds in the case of the 
DLD model as well as in the case of the BM model. In the latter case, for s < 1 or for 1 < s, the friction force is either 
infinite or zero respectively. This is due to a long range memory effect. The cutoff frequ ency u> c , by itself, does not 
prevent this feature. Finite results for the components of retarded force in the BM model ( |3.8| ) may b e obt ained only 
for bounded systems, where x(t) explores only a finite portion of the space. The calculation above (4.10) manifests 
the fact that a finite result, in case of unbounded system, may be obtained if a physical cutoff Tg is introduced. Such 
arise naturally in case of the DLD model. 

Similar picture emerges upon calculation of the effective mass. For the BM model, a finite result is obtained in case 
of an ohmic bath (s = 1). 



Am = 



00 d 



dT 



/2tt 



=e 2 



t ■ dr = —\— — for s = 1 

V 7T UJ C 

Also for 2 < s, the power law tail of a(r) decays sufficiently fast to guarantee a finite result: 

Am = — 2 s/2 T(- - 1) • nuj s - 2 for 2 < s . 

2tt 2 



(4.11) 



(4.12) 



For < s < 2 (s 7^ 1), the effective mass is infinite. In order to obtain a finite result we should turn to the DLD 
model (3.10). Here, as in the case of friction calculation, the spectral cutoff uo c is not important, and will be taken 
t o in finity. Th us the integration may be performed using the asymptotic expression for a(r). Substitution of ( |2.(| ). 
(Q and yields 



Am 



^ C0S (- s )r( s ) 

7T Z 



T 8+l dT 2 



To e 2( -t> 



T dT 



.£^cos(f s )r( s )r(i-f) 2-".Tpf 

7T Z Z 



,s-2 



[for 0<s<2] 



(4.13) 



For the special case s = 1 one obtains Am = w hich is consistent with ( 4.11 ) in the limit lo c — > oo. For s^2a 
negative infinite value is obtained, while from (4.12) a positive infinite value is obtained. In order to describe correctly 
the crossover at s = 2, one should introduce finite cutoff lu c as well as finite Tg. The calculation will not be carried 
here. 



V. PROPAGATOR FOR NON DYNAMICAL ENVIRONMENT 



In this section we shall develop a path-integral expression for the propagator of a particle that interacts either with 
static (quenched) or noisy environment. First we develop a classical expression, and then we generalize to the quantal 
regime. 
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A. Classical Derivation 



We refer to the dynamics generated by the Hamiltonian (2.7). The classical Liouville propagator, over infinitesimal 
time dr, and for definite realization of hi is: 



K{R 2 ,P 2 \R x ,P x ) = 2tv8 (P 2 -Pi) + ( 



)dr) ■ 5[ (R 2 -R x ) - —dr 



,dV dU, 
r dR + dR' 

It is more convenient to write an expression for the propagator of the Fourier transformed probability function 
(p(R,P) p(Rr)), namely, 



K{R 2) r 2 \R x ,r x ) = ( 



2iThdT ' 



exp 



If. ,R 2 -Ri. ( . ..dV dU..r 2 +r x 



Here a dummy parameter h has been inserted. Its value does not have any effect here. However, later comparison to 
the quantum mechanical version will be more transparent. For finite time, the convolved propagator may be written 
as a functional integral: 



K.(R : r\R ,r ) ■ 
Here the measure is 



R(0) Jr(0) 



VRVr exp 



im / drRr 







* dV 
dr— ■ r(r) - l 







dr— ■ r(r) 



VRVr = 



\2ttMtJ 



# segments 



.dRj,dR 2 dR x ...dr^dr 2 dr x 



(5.1) 



and the restrictions at the endpoints are -R(O) — Rq, r(0) =ro and R(t) — R, r(t)=r. It is now possible to average over 
realizations of hi, using the well known Gaussian identity 



,i J dTk(r)z(r)\ _ J drdr' (z(r)z(T')>fe(r)fe(-r') 

One obtains 

K(R,r\Bc,r Q ) = f f VRVr <•''*•-'•"•"•'• <> 
where 



_R(0)=_R o Jr(0)=r 



5 e / / [R, r] = Sfree [R, ^ ~ J dT ~QR' r ^ [classical 



and 



S N [Rr] 



2 Jo Jo 



dr x dr 2 [-<!>" (R(t 2 )-R(t x ),t 2 -t x )] r{r 2 )r{r x ) 



(5.2) 



(5.3) 



(5.4) 



(5.5) 



The free part of the effective action, for use in ( |5.3| ), is Sf ree [R, r] = m Jo drRr. It may be more convenient to write 
down the path-integral expression for K,(R, P\Rq, Pq). This expression is obtained by double Fourier transform of 
K,(R, r\R 0: r ). The result is 



Here. 



JC(R,P\Ro,P )= / VR Vr e < * s «" [ * ,r] e~P 



Sf ree [R,r] = — m I drRr [For use in (5.6)] 



S N [R,r] 



(5.6) 



(5.7) 



Note that the integration Vr is not restricted at the end-points, whereas the integration VR is restricted at the 
end-points both in R and in R. The restriction on R at the endpoints is implicit, through the dispersion relation 
it = P/m. 
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B. Qauntal Derivation 



A similar e xpre ssion may be obtained for the quantal propagator. Again, we refer to the dynamics generated by the 
Hamiltonian (2.7). The environment may be either "noisy" or "quenched", where the latter case constitute formally a 
special case of the former. The expression that will be obtained is a generalization of a result that has been obtained 
m Ref. |1J for white noise potential. 

The Feynman path-integral expression for the propagator of the quantal wave-function is, 

U(x\x )= T Vx e *tfiMt™*-v(x)-u(x,t)) (5 8) 

Jx(0)=x 

The path integral expression for the propagator of the density probability function constitutes summation T>x'T>x" over 
the pairs of paths x'(t) and x"(t). Alternatively, we may use also the coordinates R = (x'+x")/2 and r = (x"—x'), 
thus the summation will be VRDr, namely 

K(R,r\Ro,r ) = f VR f Vr e i * 5 "' [ *' r] - i * /o <M"(*V)-"(*V)) 
where 

S eff [R,r] = S free [R,r]- [ dr{V{x", r)-V{x', r)) (5.9) 

Jo 

It is important to notice that the quantal definition of the measure is ide ntica l with the classical one ( pM| ). In order 
to perform the average over realizations of U using the Gaussian identity (5.2), one may write the last expression as 



lC(R,r\Ro,r ) = f f VRVr e 

J Rn J rn 



i±S eff \R,r]-ii f Q dr f^dz(8(z-x")~8(z-x'))U(z,T) 

I R a J r a 

One easily find that the final result may be cast to the form of equation (|5.3|) or (|5.6|) with 



SN[R,r} = - / dndn {(p(4-x'l,T 2 -r 1 ) + 4>(x' 2 -x' 1 ,T2-n)-2cj ) (x^-x' 1 ,r2-T 1 )) (5.10) 

1 JQ JO 

where Xi IS cL short notation for x(j~i). 



We are now in position to compare the classical propagator ((p.q) with (5.4 )-(5.5)), with the quantal one ((5.6) with 
( 5.4 )-( 5.5 ) ) . In the latter case Ti is, in general, no longer a "dummy variable". The exception being the case where 
the actions are quadratic in the path variables, which is the case with BM model provided V(x,t) is quadratic. The 
"classical feature" may be characterized as arising from invariance under the scaling transformation of the auxiliary 
integration- variable r(r). In the quantal regime the replacement Ti — > Xh cannot be compensated by the scaling 
r — > Xr. Note however that the limit h — > is equivalent to taking leading behavior of the actions in the limit r — > 0. 



VI. PROPAGATOR FOR DYNAMICAL ENVIRONMENT 



A. Feynman Vernon Formulation 

Here we follow closely the notations in Ref. |gl| . The path- integral expression for the reduced propagator of the 
probability density function is of the general form ( |5.t| ), with 

S eff [R,r] = S free [R,r}~ f dT(V(x", t)-V(x', t)) + S F [x", x'] (6.1) 

JO 

The expressions for the reduced-action- functionals Sf[x' , x"\ and Sm[x' , x"\ in case of BM model, are given in equations 
(2.13)-(2.14) of the latter reference. In the BM model the interaction is via the dynamical variable x, while, in DLD 
model, the interaction is via u a (x). Thus, in the expressions for the friction functional Sf[x',x"] and for the noise 
functional Sn[x',x"] one should make the replacements x' — > u a (x') , x" — > u a (x") , and sum over a. Thus one 
obtains 

S F [x',x"] = den 1 gity x J J dridr 2 2&(t 2 -ti) ^ '^(u a (x2)-u a (x'2))(u a (x")+u a (x' 1 )) 
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and 



1 1 

Sn[x',x"] = dengity x - J J dTidr 2 0(t 2 -ti) ^ (tt a (a: 2 ' ) —u a (x 2 )) (m q (a//) -ti a (a^ ) ) 



Above <S(t) = a(r)0(r). These expressions may be simplified using (2.12). The results are 



S F [x',x"] = / / dridT 2 2a(r 2 -ri) • \{w{x' 2 '-x'l) -w(x' 2 -x' 1 ) + w{x' 2 -x' 1 ) -w(x' 2 -x'{)] 
Jo Jo 2 



(6.2) 



and 

S N [x',x"] = ^ / / dT!dr 2 0(t 2 -ti) [ty(4'-a#) + tu(4-ati) - 2w(4'-a/ 1 )] (6.3) 



2 



In the next paragraph we discuss further the physical significance of the functional Sf\x ' 



x 



For the BM model it has been noted before that in practice one could substitute (2.5). Thus 

S F [x',x"}= f f d n dT 2 2a{r 2 -n) r(r 2 ) ■ (R( n ) - R(t 2 )) . 
Jo Jo 

However, the formally "correct" expression is somewhat different (Ref. |2l[] equation (2.13)). Namely, the integrand 
in the above equation is with r(r 2 )i?(ri), rather than r(T 2 )(R(T±)— i?(r 2 )). One may say that the BM reduced-action 
functional includes an additional term. It is not difficult to demonstrate that the latter can be split into two terms that 
corresponds exactly to J-Apotentiai and to T SW it c hing discussed in Sec.III(B). The ^potential term may be absorbed 
in the definition of V(x), while the switching term, which is — Afc r(0)i?(0), may be factored out of the path integral 
expression (see Ref. |2lJ equation (2.34)). It has the effect of operating on the initial probability function with an 
impulse that acts on the particle if it starts its trajectory in a point x ^ 0. This term originates due to the fact that 
the initial preparation is such that the bath is in thermal equilibrium provided x = 0. Both the "switching" term and 
the additional effective potential, are absent in DLD model. This obvious results stems from the assumed inherent 



homogeneity of the environment. We turn now to the general expression (6.2). Again, as in Sec.III(B), it is convenient 
to express a(r) as a sum of its symmetric, and antisymmetric continuations. Here we focus on the resultant friction 
functional, which is 

Sf\ friction = d,Tidr 2 a(r 2 -Ti) w(x" (t 2 )-x' (n j) (6.4) 

Jo Jo 

For Ohmic bath one obtains 

S F = V [ <Jt w(t(t)) R(t) -> -n [ dr r(r) R(t) (6.5) 
Jo Jo 

where in the last stage we indicated the classical limit. Note that the classical expression for S F can be easily derived. 
For this one should include the friction force in the derivation of Sec.V(A). 



B. Derivation for Generic Bath 



A totally different derivation of the path integral expression for the propagator is possible in the general case of weak 
coupling to thermal, possibly chaotic bath. This derivation, in case of linear coupling, has been introduce already by 
FV Ref. @. The case of general, nonlinear coupling, has been considered by Mohring and Smilansky Ref. |^|. Here 
we shall take a step further, and demonstrate that under "normal" circumstances, it will reduce to an ohmic DLD 
model. We consider a bath Hamiltonian of the general form: 

Hbath + ftint = \n)E n (n\ + }^ \rn)U mn {x){n\ , (6.6) 

n nm 

where \n) and E n are the eigenstates and eigenenergies of the bath Hamiltonian in the absence of the c oupli ng. 
This Hamiltonian depends on x, the system variable, as a parameter. However, once the full Hamiltonian ( 2.1 0| ) is 



considered, x becomes a dynamical variable. The so called Influence Functional, in our notations is (by definition): 
-" s '' + ' SnI : : U mn [x"(r)] U^ n [x'(T)] , (6.7) 



e 
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Units with h — 1 are used here. U[x(t)] is the evolution operator of the bath, in the presence of the "driving 
force" x(t). The bath is assumed to be in canonical thermal equilibrium. The probability of the n-th eigenstate is 
p n oc e^ l3Erl . Using leading order perturbation theory one obtains 



U mn [x(T)} » -i drW™(i(r)) e ,(£ '"- £ " )T for m + n . (6.8) 



Similar expression holds for m=n (see Ref. 0). Substitution into (6.7) yields: 

(iS F -S N ) = J^Pnlf [ dT 2dn W mn (4')W* n (xi) e^-E^-r,) 

Lio Jo 

dr 2 dr x (u mn {x'i)U mn {x'{) e-Wn-E^-n) ^ U ; nn (x' 2 )U: nn (x'i) e' 



(E m -E n ){T 2 - T1 



(6.9) 



Now we take a further assumption which will reduce the resultant expression for Sf and Sn to the form of (6.2) and 



(S.3) respectively. The matrix elements U mn are assumed to be real, while their dependence on x is assumed to be 



characterized by the function 

n ^ (Pn -Pm) U mn (x 2 )U mn (xi) S(uj - (E m -E n )) = w(x 2 ~X\ ) ■ J(u) . (6.10) 



The reduction of fl6.9|) to ( p.2[) and ( |6.lj| ) is easily verified via 
algebraic manipulation, using 9E nm Pn{»}]= E„ m (??™-Pm){..} and ME nm p4--}]= E„™ coth(±/3w)(p„-p m ){..} 
where coth(|/3cj)=(p n +p m )/(p n — p m ). Thus, for general nonlinear coupling, DLD model constitutes an equivalent 
representation for the bath, as far as the reduced dynamics of the system is concerned. For the particular case of the 
weak linear coupling, there is further reduction to the BM model, with 

J(uj) = 7T^ (p n -p m ) \Q mn \ 2 S(uJ - (E m -E n )) (6-H) 
n,m 

In the above expression Q denotes the collective bath degree of freedom via which the interaction takes place, namely, 
in ( |6.6| ) one should substitute U{x) = Q ■ x. The bath variable Q may be some complicated nonlinear combination of 
many elementary bath degrees of freedom. For both, the DLD model and the BM model, the expression for J(ui) , 
may be cast into the form 

m - (l-e-**) a(u, f K(u) (6.12) 

where Z{0) is the partition function, cr(u>) is the standard deviation of the off-diagonal matrix elements (a; being 
the offset), and K (u>) is the density-density autocorrelation function of the spectrum {E n }, appropriately averaged 
over the relevant energy scale. The latter is related, for small energy differences, to the level spacing distribution. 
However, for any practical use, one should ignore the effect of level spacing statistics on J {to), since it corresponds to 
non-physically very long times. Thus, the generic behavior of J{uj), for physically-relevant small u, is expected to be 
linear. This leads to the conclusion that under "normal" circumstances the ohmic DLD model is good representation 
for the dissipation process. This conclusion does not hold if strong coupling to a chaotic bath is considered. In the 
latter case Zener transitions may dominate the dissipation process |],|l(j. We do not know whether the Influence 



Functional (6/7) for that case can be reduced to a form that resembles that of the DLD model. 



C. Loss of Quantal Interference 

The suppression of quantal interference is an important issue in both semiclassical and mesoscopic physics. It 
may arise that the quantum mechanical propagator may be expressed as a sum of probabilities to go either via one 
classical trajectory x a (r) or via a different classical trajectory Xb(r), plus an interference term. The expression for the 
Influence Functional may be used in order to compute the suppression of the interference due to the interaction with 
the environment. The interference term is multiplied by a "dephasing" factor (e tlp ), where we follow a notation due to 



Stern, Aharonov and Imry |18|. In our notations, the dephasing factor is identified with exp(— S , Ar[a; (r), Xb(r)]). We 
defer further discussion of interference within the framework of FV formalism to the last paragraph of this subsection. 
It is enlightening to consider the case of white noise, with 4>{t 2 — t\) = v5(t 2 —ti). For the BM model one obtains 



{e iv ) = exp 



i v r l 



^ Jo 



(x a (T) - x b (r)) 2 dr 



(6.13) 
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Thus, interference is suppressed more effectively if the two interfering paths are better separated. A totally different 
result is obtained in case of the DLD model. Here we assume that the two interfering paths are well separated with 
respect to the microscopic scale £, namely £ <C \x a (T)— Xb(r)\ most of the time. It follows that 



3 iV ) = exp 



2 



■ t 



(6.14) 



Here, the interference decays exponentially in time, and the actual spatial separation of the paths play no role. Due 
to the disorder, dephasing events are as effective for small separations as for large separations. 

The dephasing of interference in metals due to electromagnetic fluctuations has been discussed by APtshuler, Aronov 
and Khmelnitskii in Ref. Jl^]. Their results have been re-derived by Stern, Aharonov and Imry [j^). A somewhat 
simplified derivation is reconstructed in App.A. The strong dimensionality dep enden ce of the deph asing process has 
been emphasized. This dependence is due to the separation-dependence as in ( 6.13| ), and see also (|A.(j ). In case of 



the DLD model, the local nature of the dephasing process will eliminate this feature. 

In order to understand how interference arise from the FV path integral expression ( |5.6| ) , it is convenient to rewrite 
it in the following form 

i-R.p 

JC{R, P\R , P ) = VR JC[R] , (6.15) 

•ARo.Po 

where IC[R] is a real functional, which is defined by the expression: 

K[R] = f Vr e z i s "ff^ e^ 5 "^ . (6.16) 

J unrestricted 



If the path integral expression for the evolution operator (5.8) is dominated by a single classical path x a (r), then 



the T>R integration in (6.15) will be dominated by R(t) = x a (r). Obviously, in order to obtain a non - vani shing 
result, the endpoint conditions should be compatible. Turning to the computation of /C[i?(r)=x a (r)] via J6.16| ), one 
observes that the T>r integration is dominated by the trivial trajectory rfl(r)=0. We use the subscript R in order 
to suggest that in general this trajectory should be independent. Indeed, this is the case if two classical trajectories 
x a (r) and Xb(r) dominate. One should consider then the "interference path" R(T)=(x a (T)+Xb(T)) /2, for which the 
T>r integration is dominated by the non-trivial paths t\r(t)= ± (x a (j)— Xb(r)). The existence of non-trivial path rn 
is the fingerprint of interference phenomena. A classical trajectory R(t) for which r^=0 will not be damped, since 
SN[R,rn] = then. In contrast, an interference path, for which t\r ^ 0, is damped, since < Sn[R, tr] in general. 
However, in Sec. IX, where localization effect is discussed, we shall encounter a vast family of interference trajectories 
that are not damped by the noise functional. In the latter case, the interference paths are found to be dominant in 
the computation of the propagator. 



Another issue that deserves attention is the interplay of friction and interference. By inspection of (6.4) it is clear 
that for disordered environment, quantal interference is unaffected by friction. This is true as long as x a and Xb are 
well separated in space. In Sec.VII(B) we shall encounter a related quantal manifestation of this observation. 



VII. SPREADING AND DIFFUSION 



The main results of the two last sections are the path integral e xpr ession (|5 
with the app ropr iate action functionals Sf ree ([sTrj), Sf\ friction (5.4), and Sjv 
presented in fl5.5|). For ohmic friction, both the quantal version and its classical limit (6.5) will be further considered 



) for the propagator K{R, P\Rq, Pq), 
3|). The classical limit of Sn is 



and compared. Friction in case of non-ohmic bath has been discussed in Sec. IV, and its quantal analog will not be 
considered further. 

In order to get preliminary insight into the path-integral expression, co nsid er first the case of free part icle in "white" 
noisy environment. Namely, <p(r— r') = v8{t—t'), with v = 2rjkBT as in 4/7. In the classical case (5.5), one obtains 



1 f* 
SN[r] = r J Q 



r(T) 2 dr 



(7.1) 



independent of the spatial auto-correlation function w(r). The observation that spatial correlations are of no impor- 
tance, as long as the noise is uncorrelated in time, is trivial from classical point of view. In the quantum mechanical 
case, the corresponding expression is 



SjvH = v / (w(0)— w(r(r)) dr 



(7.2) 
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In contrary to classical intuition, spatial correlations may be of importance. However, for the BM model, ( |2.5| ), one 
recovers the classical result. 

In the path integral expression ( |5.6| ), one may perform the integration Dr = ...dr^d^dri . In the absence of Sn, 
each integration over dr n results in delta function of velocities, namely <5(i?(r n+ i) — R(r n )). In the presence of Sn, 
the integration dr n is weighted, and as a result, each S function is smeared. The propagator constitutes a convolution 
of these smeared S functions. In particular, both in the classical case and in the BM model, each smeared 5 function 
is a Gaussian. It is obvious, that both in the classical and in the quantal case, Sn leads to stochastic-like spreading. 
In what follow we want to estimate this spreading. 



A. Non Disordered Environment 



We turn now to estimate the spreading in the absence of disorder, which is the standard BM model. The noise 
functional is quadratic in the path-variable r, and is independent of the path-variable R, namely 



1 



o jo 



SN[R,r] = -l I dndn <P(r 2 -ri) r(r)r(<0 (7.3) 



Here, an exact treatment is available Q . One may expand both Sf ree and Sn around the so-called classical paths, that 
are determined by the variation 5Sf ree — 0, with the constraints 5R—Sr=0 at the endpoints. The Gaussian integration 
is performed exactly. General expressions may be found in Ref. pi]] . The phase space propagator JC{R, P\Rq, Pq) is 
obtained by double Fourier transform of /C(i?, r\Ro, ro), and obviously results in Gaussian function. In particular, one 
is interested in the spatial spreading. Setting Pq = and integrating over the final momentum P, one obtains 



K(R\Ro) = - — exp 

Z7T^/ G spatial 



1 / R-R, 



2 \ O spatial , 

A general expression for the spatial spreading, that applies to any </>(r), may be obtained p] 



(7.4) 



o- S patiai = — \ [ [ 4>{t-t') ^(t) r c t(r') drdr' , (7.5) 
m \ Jo Jo 

where r c i solves the linearized classical equation of motion (mr + rjr = 0) with initial conditions r(0) = and r(0) = 1. 
Note that ohmic friction is considered, for which the BM model is well defined. For the simplest case of white noise 

without friction one obtains a spa tiai = \l \ v ~$p- Friction leads to damping and diffusion. In the latter case, considering 
long time e - ™' -C 1, one may disregard a short transient and substitute r c ; = mjr\. Consequently cr spatia i = ^ jpt. 
However, at low temperatures 4>{t) = — while <fi(T)dT = 0. Thus 



la l 




special case - suppressed diffusion , (7.6) 



diffusion is suppressed due to the negative autocorrelations of the noise. This effect is classical in nature. Intentionally 
we did not use the explicit expression for the constant, namely C = hrj. The presence of h~b a th in the formula, rather 
than Tisystem, without explicit subscript, may miss-lead the reader. The particle can be treated as a classical object, 
within the framework of e.g. Langevin equation, and still the suppression of diffusion will occur. 



B. Disordered Environment, White Noise 



Both t he classic al li mit and the BM model generate "classical" dissipation effects. We turn now to the DLD 



model, ((7.2) with ( |2.6| )). Here the situation is quite different. In order to compute the propagator, a different, more 
powerful strategy is required. We shall exploit the fact that for white noise Sn[R,t\ is still independent of R(t), 
while S e ff[R, r] is linear in R(t). It is most convenient to Fourier transform /C(i?, r\Ro, ro) to lC(p, r\po, ro). The path 
integral expression may be written for this representation as follows ( we suppress from now on printing h~, but shall 
restore it later): 

K,(P, r\po,r ) = e ir l(R™'{r)-R w'(r )) f'^. f VR & -i f Q dr( m f+riw"(r)r)R & -v J*[w(0)-w(r(r))]dT 
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The T>R integration can be performed now, yielding a delta function 5(mr+r)w" (r)f) at each point along the trajectory. 
The trajectory for which the integrand does not vanish will be denoted by r(ro, r, t), where < r < t. This trajectory 
is required to satisfy the endpoints conditions r(0)= tq and r(t) = r. Hence, the result of the path-integration is 



K(p,r\pa,r ) 



— c irt{Rw' (r)-R w' (r )) 



6(p- 



■I j-W r/ -I n\\ — v V [w{0)—w{r(r a ,r,T)]dT 

mr{r ,r,t)) 0{p — mr(r , r, 0)) ■ e Jo 



The Inverse Fourier transform yields 

K,{R r\R r ) = e i [(. m Hro,r,t)+r l w' (r))R-(mf(.ro,r,0)+nw' (r ))Ro] e ~ v j Q [w{Q)-w{r{r Q ,r,T)}dr 

To simplify the latter expression we note that mr+nw(r) is a constant of the motion. 

The propagator IC(R, P\Rq, Pq) is the Fourier transform in the variables r and r$. In order to get insight we restrict 
ourselves to the reduced kernel K,(R\Rq), namely 



K(R\Ro) = 



dr ^[(^(^fi^+vw' {r ))(R- R )} e~ u Jg( w (°)~ w ( r ( r o^. T )) dT 



(7.7) 



As before we distinguish the case of frictionless propagation, for which r(ro, 0, r) = (t—r) /t-r , from the case of damped 
particle with e~ ™* < C 1. In the latter case the trajectory is modified for |ro|<£, where r(ro,0,r) « r$. Consequently 
the "phase" in (7.7) is S e f f(ro)=—m/t-(R—Ro)-r and S e ff(ro)=—r]-(R—Ro)-ro for the two corresponding types of 



trajectories. As for the noise argument in fl7.7| ), it is Sn^o) = vi 2 t for £ -C |r |, and Sn(tq) » or ^n(i~o) ~ 5^ r o 

for |ro|<^, depending whether friction is absent or present respectively. The separation of scales both in S e ff(ro) and 
in Sn(to) enables splitting the integral in a convenient way, namely 



drn 



W 



(?) *» 



1 - W 



(?) 



dr Q 



where W(x) is a s mooth, symmetric cutoff function that equals « 1 for \x\ < 1 and equals 
integration in (7.7) is performed, and the following expression is obtained for the propagator. 



for 1 < |a;|. The 



tC(R\Ro) 



W I - 



R 



\h/r)l 



* tCd(R- Rd(t)) + 



S(R-Ro(t)) - W 



R-Rq(t) 
fit I ml 



(7. 



The symbol * stands for convolution, and results in smearing of the classical propagator on scale Ti/ (t)£). For frictionless 
propagation one should use the replacement h/(rj£) — > M/(m£). We have used the notation Ro(t)=R c i(t)=Ro. The 
above result holds also if Pq ^ 0, in this case Ro(t) will propagate as if friction is absent, while R c e(t) will propagate 
as in the classical limit. Obviously, if fr ictio n is indeed absent, then again Ro(t) and R c e(t) will coincide. The kernel 
K c t{R\Ro) denotes the classical result, (7.4). 

The expression for the quantal propagator demonstrates that a piece of the wavepacket is frozen due to the disorder. 
This is a non trivial quantal effect that indeed can be entitled "Quantum Dissipation" . We emphasize again that 
such quantal effect is absent in the BM model. However, the expression for the propagator also demonstrates that 
the "quantal correction" goes to zero exponentially in time, as in the case of interference discussed in Sec.VI(C). 



VIII. CLASSICAL NON-MARKOVIAN EFFECTS 



The disadvantage of the treatments that have been presented in the preceding section, is the difficulty to extend them 
to the general case of disordered environment, namely if the disorder (the noise) is correlated in time. We therefore 
turn to a somewhat more heuristic approach, tha t will enable approximated treatment. F or t he computation of 
IC(R, P\Ro, Po) we shall use the classical limit (5.5) of Sn[R, r]. The "quantal correction" in (7.8), is not considered 
again in the present sect ion. 

In the classical limit ( 3.15 ) constitutes a formal solution of Langevin Equation. The real functional JC[R] has a 
simple probabilistic interpretation. For ohmic friction it takes the form 



K[R\ = / Vr 



' Io dT Ah ( t ) r ( T ) e ~5 Jg Jg dTdT ' 4>b.{t,t') t{t) r(r') 



(8.1) 



where 



-w"(R(t)-R(t')) ■ 0(r-r') 



•2) 
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and 



A fl (r) = mR(r) + V R(T) (8.3) 
Formally, the unrestricted Dr integration may be performed exactly, yielding the result 

K[R] = v/det[$ K ] • exp (^\ j J q drdr' $ fl (r,r') A r (t) A*(t')) (8.4) 

where $^ is the reciprocal of 4>r. In order to compute IC(R, P\Rq, Pq) one should identify the most contributing paths, 
for which JC[R] is maximal. In subsection A, where we discuss short-time correlated colored noise, we assume that one 
optimal path dominates the computation. In subsection B, where we discuss static or almost-static noisy potential, 
we shall identify a whole family of optimal paths. In the next section, where quantum localization is discussed, we 
shall use a similar strategy, and a family of interference paths will be identified. 



A. Normal, Dissipative Diffusion 



In this section we shall analyze the diffusive behavior which is encountered in the absence of disorder. Our as- 
sumption will be that the T>R integration is dominated by one smooth "optimal path" . We shall substantiate this 
assumption by demonstrating consistency with the exact result that has been presented in Sec.VII(A). Furthermore, 
it will be argued that the "optimal path" for short-time correlated noise is the same as for white noise. 

In order to find the path R (t) that maximize /C[i?], we consider first the case of white noise, where 4>r{ t ~ t ') = 
vr8(t— t'). Hence 3>r(t—t') = 6(t— t'), with vr = v (the subscript R is reserved for later use). As in Sec.VII(A) 
we focus the attention on the comp utation of the reduced propagator )C(R\Rq). Formally the path integral expression 
for JC(R\Rq) is identical with ( |6 .15 ), except for the restriction at the endpoints. For K,(R\Rq), the relaxed constraints 



are R(0)=Rq, i?(0)=0, and R(t)—R. Denoting R — v, the variational equation for i? Q (r), including Lagrange 
multiplier, is 



t 

2 



S / [(mi + nv) + const ■ v]dr = 
Jo 

one obtains 



mv — rj 2 v = const 
(mv + r)v)t = 
u(0) = 

/ v{r)dT = (R-Rq) 
Jo 



The last two equations are the constraints. The solution for damped propagation (e r»* <C 1) is easily found. For 
sake of comparison also the solution for frictionless propagation (n — 0) is displayed. 

R (t) w — ^1 — ie~m( t ^ r ) — (i — ie-^^c - ™ 1 "^ for damped propagation 
^ 

R (t) — — — ^ — r(2t — t) for frictionless propagation (8-5) 

In the first formul a a constant prefactor that equals w 1 has been dropped, since we assume here e - ™* -C 1 . The 
optimal paths of (3.5) are illustrated in Fig.||. Computation of Ajj(t) for the optimal path is str aigh t forward. For 



damped propagation the computation is trivial since Ar w t\ t is constant. Substitution into (8^) yields 



K[Ro] = const -exp ( -- ■ ( ^ - ) , (8.6) 
V 2 TP V R 1 J 

which is in consistency with the exact result of Sec.VII(A). It is easily verified that also for frictionless propagation, 
consistency with the exact result is maintained. 

For short-time correlated colored noise (1 < s) it is natural to replace <j>(r— r') by v6(t—t'), with the effective 
white noise intensity v = (f)(ui=0). As long as < v (finite temperatures), the long time behavior is diffusive, and 
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FIG. 3. Illustration of the optimal path (8.5) for either damped propagation (solid line), or frictionless propagation (dashed line). The 
horizontal axis represents spatial position, while the vertical axis is for the time. 



consistency with ( J7.5| ) is easily verified. Still, a more elaborated argument is required in order to substantiate the 
"white noise approximation" . This argument will be discussed now. 

By inspection of ( |8.4| ) it is clear that the most contributing paths, for which fC[R] is large, must satisfy |A^((jj)| 2 < 
</>(w). For white noise <j}{uj)=v=const 1 but still the "optimal path" is smooth. By "smooth" one means that A^(oj) 
is concentrated within the interval uj < 1/r^ where T v =m/ri is the relevant time scale for the system's dynamics. 
Consider now a noise auto-correlation function of the form <p(r) — v8(t) + C-G s (t). Its Fourier transform satisfies 
4>(uj) ~ v w const for lu < {v /C) 1 ^^ 1 ^. Thus, the first requirement for the "white noise approximation" to hold 
should be C/t^ 1 < v. For oj c < u> the spectral function 4>(oj) drops to zero, which implies that such high frequency 
components are not favored. So far there is consistency with our assumption that the most contributing paths are 
smooth. However, in the vicinity of lu c the spectral function 4>{uj) is peaked. It implies, that unlike the case of white 
noise, an oscillatory component with time period r c is favored. Obviously, such component arise from the strong 
accelerations that the particle experiences within short perio ds w hose duration is r c . Over these short periods the 
maximum displacement is AL = i(y / (/)(r=0)/m)r 2 . Using (^4) it is found that this amplitude is proportional to 

tJ 4 ~ s ^ 2 . For s < 4 the amplitude AL goes to zero as t c — * 0. Therefore, in this restricted regime (1 < s < 4), and in 
particular for s — 2 (low temperature ohmic noise) the "white noise approximation" should be adequate. 



B. Anomalous "Diffusion" 



Encouraged by the consistency of the heuristic approach with exact results, we turn now to analyze the diffusion 
due to short-time correlated noise in the presence of disorder. We shall use the "white noise approximation" whose 
validity has been discussed in the preceding subsection. Namely, for short range correlated no ise, the most contributing 
paths are concentrated around the same R (t) that has been found for white noise. From ( |S.5| ) it follows that R (t) 
is, up to end point transients, a free-like propagation. Consequently, the effective noise autocorrelation function is 
(I>r{t) ~ —w"((R—Ro)/t ■ t) • (j){r) and we define 

v R = J°° [-w" (^t)^(t)] dr (8.7) 

In general < vr, also in the limit of zero temperature. 

In the presence of disorder, the effective white noise intensity is (in general) a function of the endpoint conditions, 
rather than a constant. For typical noise autocorrelation function of the form 4>{t) = C ■ G s (t) one obtains 



with ti 



i-t 

(R-Ro) 
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and 



C = 



(8.8) 



The computation has been carried out by taking in (3.7) the Fou rier transform of both </>() and w"(), and performing 
du> integration rather than dr integration. Substitution into (3.6) suggests that: 



K{R\Rq) \{\R-R \<l Wc t} 
K,(R\R ) |{|Jl-.Ro|><W c t} 



const ■ exp 



const ■ exp 



4f 

2C 
V 



1 t s - 2 \R-Ro\ 3 - s 



\R-Ro\ 5 



2C LVc +Z £ 3 t 4 



(8.9) 



The "tail" of the dispersion profile is universal. It depends on cj c , but it is independent of the nature of the noise. 
In contrast, the short range profile is determined by the low frequency bath-oscillators, and thus it is sensitive to the 
exact value of s. For s = 1, (ohmic model, high temperatures), normal diffusive behavior prevails. For s = 2 (ohmic 
model, low temperatures) the diffusion freezes. The dispersion profile is Exponential rather than Gaussian, namely 



K(R\Ro) 



const ■ exp 



V 



\R-R 



2 1 f~i 

7T 7^1 ^ 



[ohmic, low temperatures] 



(8.10) 



The dispersion here, in the DLD model, is of the order C/(n 2 £) rather than \fC jr\. (BM model (7.6)). One should 
observe that on for C ~ (j]£) 2 both models predict consistently dispersion on spatial scale £. For larger noise intensity, 
the BM model is not valid, and DLD predicts always a larger dispersion, which is intuitively expected. For weak 
noise (C < (ri£) 2 ) the spatial spreading is on scale less than I. Our treatment of the DLD model is not valid on this 
microscopic scale. However, in this regime the BM model can be trusted. For 3 < s < 4 equation (3.9) implies that 
the particle is evacuates from the vicinit y of Rp. 

The "body" of the "diffusion" profile (8.9a) is determined by a low velocity paths for which r c -C Tg. It is easily 
verified that a sufficient condition for the validity of the "white noise approximation" is rg <C T^. This condition 
is satisfied in the relevant spatial range (\R—Rq\ < £u> c t), except for a relatively small interval around Rq, which is 
determined by the large ratio T^/r c . 

The "tail" of the "diffusion" profile (3.9b) is determined by a high velocity paths for which rg <C t c . Here the the 
validity argument should be modified. The spectral function <p(w) is peaked around u> = 1/rg rather than around 
lu = lu c . Thus, oscillatory component which is characterized by period Tt is favored. The maximal spatial amplitude 
of this component is AL = i(y / <^(r=0) /m)r| . It is convenient to use a special notation for the standard deviation 

\{W/mv 2 )£. 



of the disordered potential, namely W = £a/ 4>{t- 
This amplitude is required to be much less than i 



0). Using this notation the amplitude is AL 



or alternatively W <C 



or alternatively ti <C tw, with 



t~w — (/ y/2W/m. In this subsection we have limited the discussion to the case of short-time correlated noise for 
which t c is small in some sense. Indeed, if it is assumed that r c <§; tw, then the validity condition will be satisfied 
automatically. 



C. Classical Localization and Non-Dissipative Diffusion 

In this subsection we shall discuss the case of static or almost-static disordered non-dynamical environment (rj = 0) . 
The "noise" is assumed to possess long-time correlations. Specifically: 

W 2 ( Kt-t') 2 \ 

Mr,r') = ~w"{R{t)-R{t')) ■ —expU-i— ± J , (8.11) 

where W is the s tandard deviation of the disordered potential. Here r c is assumed to be large, much larger than 
t w = £/y/2W/m. 

In general the "white noise approximation" breaks down. T he p ath integral expression is no-longer dominated by 
"smooth" trajectories. It is difficult to use the explicit formula ( |S.4| ) in order to identify the family of "optimal paths" 
since $(r, r') is no longer diagonal. We therefore prefer to use heuristic considerations. No new insight is gained if 
one insists using (8.4). 

We take first the limit r c — > oo. The classical mean free path of the particle is 

/ 1 \mv 2 \ 

^hV ■ (8 - 12) 
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The corresponding time is T co m S i on = L collision/ v. After that time the probability of being backscattered is of order 
1. In one dimension this backscattering will lead to (classical) localization of the particle. The localization length is 
exponentially large for high energies. 

If t c is finite, rather than infinite, classical localization will manifest itself only if T co m s i on <C t c . Within the time 
scale t c , the particle spreads over spatial range of the order L co iiisiom while its velocity is randomized. It follows that 
K,(R\Rq) may be used as a stochastic kernel. Hence the Markovian property is recovered over time scales that are 
much larger than r c , and a diffusive behavior follows with coefficient L 2 olHsion /T c . This diffusion is non-dissipative 
"random walk" like. 



IX. QUANTUM LOCALIZATION 

As in the last subsection, we shall discuss the case of quenched disorder. However, here the quanta! analysis will 
be carried out. The v aria nce of the disordered potential is W 2 , with auto-correlation length I. We approximate the 
Gaussian correlation (|2.6|) by a delta function. Defining a = the path-integral expression for the propagator is 



t2 



fl5.6| ) with ([5J) and 

S N [x',x"} = \aW 2 I I (5(x%-x'{) + tfOa-z'i) - 25(x%-x' 1 ))dndT 2 (9.1) 
2 Jo Jo 

The essential feature of this functional, is its non-local nature. Our first step will be to get some insight into Sn- 

Let us consider segment i and segment j that belong both to the same path, either x' or x" , within the same spatial 
interval Ax. This includes the possibility i — j. The contribution to Sn is 



ASV = + -aW 2 



2 \viVj\ 

where v is the velocity x within Ax. However, if i belongs to x' , while j belongs to x" , or vise versa, then the 
contribution to Sn is 



Ax 

N 

\ViVj 



AS N = -aW 2 



One easily convince oneself that the contribution of each spatial interval Ax is non-negative. A zero value may be 
obtained if to each segment i, that belongs to x', corresponds segment j that belongs to x", with \vj\ = For 
example, a zero value for Sn may be obtained if x"(r) is shifted in time with respect to x'(r). Referring to the 
endpoints, one should identify r = t with r = 0. It follows that Sn = implies that the two paths x' and x" are 
either identical, or satisfy the constraint R(t) = R(0). If R(t) ^ -R(O), the value Sn = may be obtained only if x' 
and x" are identical. More generally, if \v t \ = \vq\, one may prove that the following inequality holds for any pair of 
smooth paths x' and x", 

e -^S N [R,r] < e -(f^mi„(|ii t -i? |,|(|r t | + |ro|)) _ (Q2) 

For particular i?(r) one may ask what is the r{j) for which Sn[R, r] is minimal. The trivial minimum, which is 
also the absolute minimum, is r(r) = 0, for which SN[R,r] = 0. However, any small perturbation on r(r) will 
make Sn[R, r] much larger. Therefore, we are tempted to assume that there may be some other, more stable (local) 
minimum ?\r(t). A non-trivial local minimum r^(r) does not exist for any R{t). However, one can prove that there 
is a large family of R-s for which such minimum exists, by actually constructing them. This is done by following the 
considerations that were presented at the beginning of this paragraph, An example for such construction is presented 
in Fig.||. The situation here should be contrasted with that encountered in Sec. VIII. In subsection A (there) we could 
have defined one optimal path R (t). Here, there is a whole family of "optimal paths" , as in the case of subsection B 
(there). However, in the present case these paths are "interference paths" rather than "classical paths". The following 
observations concerning the relevant optimal paths are important: (a) They consist of many straight segments, and 
have ZigZag character; (b) The final velocity R(t) is favored to be equal in absolute value to the initial velocity; (c) 
Turning points impose significant restrictions; (d) The non trivial minimum r^(r) is isolated; (e) The non-trivial 
minimum is relatively stable. The last point is the most difficult to observe. First it should be noted that (e) must 
be true a-priori. Else, if the trivial minimum dominates the path-integral expression, then the result would be that 
the particle has roughly the same probability to go from any initial conditions to any final conditions, irrespective 
of proximity considerations or even energy conservation. Still, a reasonable argument is required why the non trivial 
minimum is relatively stable. For this consider a straight segment i for which R(t) = Vi and tr{t) — ri. One observes 
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FIG. 4. Illustration ol an optimal path R (solid curve) for localization problem. The dotted lines are x'(t) and x"(t) that correspond 
to rn(r). The horizontal axis represents spatial position, while the vertical axis is for the time. 



that if r R is perturbed by a fluctuation of time period |j"j|/|u|, or by some higher harmony, then the contribution to 
Sn[R,t] is negligible. This is to be contrasted with the case of the trivial minimum r(r) = 0, where any fluctuation 
has high cost. 

Wc turn now to the formal extension of the procedure that has been presented in the previous section. We expand 
5/v[i?,r] around the non-trivial minimum: 



oo i>00 



S N [R,r}^ S N [R,r R } + / <W ^(r-r') • (r(T)-r R ( T )) ■ (r (/)-?•*(/)) + ... (9.3) 
Jo Jo 

Here we are not able to write an explicit expression for the highly complicated kernel 4>r(t—t'). However, we proceed, 
and write down the result of the Gaussian integration: 

K[R] = const- cos (J dr A R {T)r R { T )^ e - SNlR ' rit] ■ 

■ exp f-i J J drdr 1 $ r (t - r') A r (t) A fl (r')) (9.4) 
Now, we should perform the VR i ntegration. This integration will be dominated by the family of "optimal paths" . 



Note that the cosine term in (9.4) equals unity since for r R the trajectories x' 
with respect to the other, hence S e ff = 0. Within the family of optimal paths, not all have the same contribution. 
One should expand around those that have the largest contribution. For these paths Sn[R, r] = ^-^-\R — Rq\. Here 

we consider endpoints conditions (i? n , mvg) at t = and (i?, ±mi>o) at time t. The time t is assumed to be sufficiently 
large to guaranty steady state distribution. It follows that 

K{R,±mv Q \Ro,mv Q ) ~ exp (-2 ^ '~ R f \ (9.5) 

where 

avo) = 2-^ . (9.6) 
aW z 

For convenience fi has been restored in the latter formula. In the application of the inequality (^^) a sub-family 
of "optimal paths" has been ignored, for which \r R {endpoints)\ < \R—Rq\. It is justified provided this sub-family 
constitutes a zero fraction of the whole family. For this we should assume sufficiently long time (t) , for which a steady 
state distribution is attained. 

Both the validity and the applicability of the inequality ( |9.2| ) demonstrates the vulnerability of our localization 
argument. It is important to consider circumstances in which either of these conditions is not satisfied. The inequality 
(|9.2[) will not be valid if the noise is not infinit ely correlated in time, i.e. the disordered potential is not static when 



viewed on large time scales. The inequality (9.2) will not be applicable if additional white noise is added to the 
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Hamiltonian. In the latter case, large tr(t) will be suppressed by the corresponding additional term in the noise 
functional. Consequently, the sub-family of paths for which \rn(endpoints)\ < \R—Rq\ will not constitute zero fraction 
of the whole f amil y. 

Expression (3.6) for the localization length agrees with the well know result for ID localization of "free" particle 



(Thouless Ref. |Tj). Note that Thoulcss uses some scaled units, resulting in the expression £ = 8mE/h 2 . We shall 



show now that the above expression is somewhat more g ene ral, and applies also to cases where the dispersion law is 
different from v = p/m. For this one should replace in (5.8|) the kinetic term 1 ■ 
resulting in 



by some general function T(x), 



dT{T(x")-T(x')) 



dr p(R) 



(9.7) 



where p(v) = T'(v) is the dispersion law. The subsequent formalism is easily generalized. The expression (9.6) for £ is 
unaffected. Actually, the general expression for £ could have been guessed. We use naively the Born approximation, 
calculate the mean free path and rely on the fact that £ is twice the mean free path (Thouless Ref. [14]]). By the 
golden rule, the probability of being backscattered is 



Prob = 2n\(-p\U\p)\i-^ 
where L is the length of the available space, and E is the kinetic energy. The matrix element is 



1 



1 



\{-p\U\pW = j- 2 -|FT[W]| 2 = -aW 2 



(9.8) 



(9.9) 



where in the last equality FT[] denotes Fourier Transform, and hi is assumed to be uncorrelated in space 
("white spatial noise"). The result (9.6) is easily recovered. One should use £ — 2£, and vq — dE/dp, 
while I is evaluated via the relation Prob = vq/£. In order to further demonstrate the generality of (p.Gj), 
let us consider Anderson tight binding model. The spacing will be denoted by a. The Hamiltonian is 
H = + T(\n)(n—1\ + |n)(n+l|)) where V n is uniformly distributed in [—Wo, Wo]- The transition am- 

plitude is T . The Kinetic energy in the center of the band is E — 2Tap, where p is the momentum. Hence vq = 2Ta 



there. The dispersion of the on-site energies is W 



function is aW 2 . Formula (9J3) suggests that £ 
(66) there), including the prefactor. 



n-2 



2Wq. The "area" under the triangular autocorrelation 



2Aa(T/Wo) ■ This results agrees with that of Ref. |13| (equation 



X. SUMMARY AND CONCLUSIONS 



A unified treatment of Diffusion Localization and Dissipation (DLD) has been presented in this work. All these 
phenomena may be derived from the general path integral expression ( |5.6[ ), 

pR.P p 

K(R,P\Ro,P )= VR Vr e l i s ^f^ e~l^ SN[R - r] , (10.1) 

JRo,Po J 

upon inclusion of the appropriate functionals S e ff and 5jv- General expressions for these functionals are available 
and various limits may be considered: (a) Quantal versus classical expressions; (b) Disordered versus non-disordered 
environment; (c) Dissipative versus non-dissipative environment; (d) Quenched versus noisy environment. In the 
classical limit the DLD model constitutes a formal solution of Langevin equation. The classical limit may be obtained 
by linearization of the quantal S e ff[R,r] with respect to r, while expanding Sn[R,t] to be quadratic in this path- 
variable. The disorder or its absence depends on the choice of the spatial auto-correlation function w(r). The 
dissipation is turned on if the friction functional Sj? is included in S e ff. The nature of the noise, whether it is 
"quenched" , "colored" or "white" , is determined by the noise kernel (j)(r) . In the latter case any combination may be 
considered as well (see further discussion at the end of this section). 

The classical BM model is well defined in terms of an appropriate Langevin equation only in the case of an ohmic 
bath. This is not the case with the classical DLD model. The classical dynamics in the latter case is well defined in 
terms of appropriate Langevin equation also for non-ohmic bath. Explicit expressions for the friction force, and for 
the effective mass have been derived. Another nice feature of the DLD model is the absence of "switching impulse" . 

Once the noise auto-correlation function 4>(t) is specified, the BM model is in-distinguishable from its classical limit. 
As long as the external potential V(x) is quadratic (at most), the quantal propagator is identical with the classical 
one, and Langevin equation can be used in order to describe correctly the time evolution of Wigner function. All the 
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quantal effects that are associated with the standard Zwanzig-Caldeira-Leggett BM-model are (formally) reproduced 
by solving the classical Langevin equation with an appropriate noise term. The DLD model is different. The non- 
stochastic, genuine quantal features of the DLD model have been discussed. These features constitute a manifestation 
of either the disorder or the chaotic nature of the bath. 

For either noisy or ohmic environment both the BM model and the DLD model leads to either spreading or diffusion. 
For the BM model, the spreading and the diffusion profiles are described by a Gaussian distribution. For the DLD 
model one should include a quantal correction. The disorder freezes a piece of the wavepacket, letting it to propagate 
as if it were a free particle. Both, the "quantal correction" to the propagator, as well as any other interference 
phenomena, die out exponentially in time. This exponential decay, due to dephasing, is independent of geometry. It 
should be contrasted with the results for loss of interference in the presence of BM-like environment p7|Jlq 1 . Another 
important observation is that for disordered environment, quantal interference is unaffected by friction. 

On the classical level it is fascinating to analyze the diffusion profile in the presence of disorder. For the low 
temperature ohmic BM-model, it is found that diffusion is suppressed, though its Gaussian profile is maintained. The 
DLD model, in the same circumstances, leads to an exponential profile that does not change with time. This new 
effect is due to the interplay of the temporal (negative) autocorrelations of the noise with the spatial disorder. Even 
more fascinating "diffusion" profiles are found for other types of noise autocorrelations. 

Quenched disorder in one dimension leads to classical localization, as well as to quantal localization. The former 
is characterized by a localization length which is exponentially large at high energies. Quantal localization on the 
other hand is dominated by interference phenomenon. We have identified the interference paths within the framework 
of FV formalism, and demonstrated how the well known exponential profile emerges. The localization length in the 
quantal case is proportional to the square of the velocity. The applicability of this result to other dispersion relations 
has been pointed out. 

It is obvious from the derivation, that localization cannot be argued if the noise is not strictly static (e.g. slowly 
modulated). Diffusive behavior is recovered if the noise possess long but finite auto-correlation time. Also the case of 
white noise "on top" of the static disorder will evidently lead to diffusion jl9| [p3l. In this latter case, which has not 
been considered in this paper, non perturbative effects may manifest themselves|20|] as in the case of the "Quantum 
Kicked Particle" (2j|. The DLD model should also account for diffusion in the presence of both quenched disorder, 
noisy potential and friction (all together). There should be a way to derive systematically well known heuristic results 
that corresponds to hopping, or variable range hopping p"i"[ . 

Classical non-dissipative "random walk" diffusion has been discussed in the restricted case of long-time correlated 
noise. It is essential to generalize the DLD model to more than one dimension in order to account for this phenomenon 
in the case of strictly static disorder. Auto-correlations of the disordered potential, both in time and in space, should 
be considered, in order to generate the dynamics which is described by Boltzmann transport equation. In the limit 
of quenched disorder, localization effect should be encountered. 
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A. INTERACTION WITH EXTERNAL BATH THAT CONSISTS OF EXTENDED FIELD MODES 



This appendix illustrates how the unified formalism that has been presented in this paper should be modified 
in order to deal with external bath that consists of extended field modes. Unlike the case of the BM-model, this 
modification is not as immediate as the mere substitution of appropriate w(r). In the present derivation the bath 
Hamiltonian is not specified, and also the assumption u a — u(x—x a ) is altered. 

For the standard derivation of the DLD model it has been assumed that the interaction Hamiltonian is ( |2.8D with 
u a = u(x—x a ), leading to the factorized noise auto-correlation function ( [2.41 ). In the more general case the classical 



deriva tion as well as the quantal derivations lead to a stochastic force that satisfies (2.3) or to the noise functional 



(5.10) respectively, with noise auto-correlation function which is 



(x-x',t-t') = ^<j) a (t-t') u a (x)u a (x r ) . (A.l) 



Here 4> a (t—t')=(Q a (t)Q a (t / )) eq , and the interaction Hamiltonian (2.8) is used with c a =l without loss of generality. 



The fields mode are still assumed to be decoupled, but the bath Hamiltonian is not specified. Instead, one relys on 
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Fluctuation-Dissipation theorem as in Ref. |1| , writing the general expression: 

f°° duj 1 
$a(t-tf) = / —J a {ui) hcoth(-hf3uj) cos[u(t-t')] . 

. n TT 2 



(A.2) 



Note that if the field modes are simple harmonic oscillators, then one should substitute J a (uj)=^^ L ^-5(uj—uj a ). 
Alternatively, one may speculate J a (uj) using known response characteristics of the bath. 



To make further progress the interaction Hamiltonian (2.8) should be further specified. Using "standing waves" 
decomposition it is assumed to be 



'Hint = ^2 (Qic cos(fc Q • x) + Q 2a sin(k a ■ x) 



(A.3) 



Substitution of the appropriate u a into ( Al), and converting the summation to an integral over directions and over 
k, the result can be cast into the following form 



(x-x',t-t') 



dk 



cj)(k,t-t') (cos(fcft- (x-x')))q 



(A.4) 



The average over directions can be performed leading to cos(k\x—x'\) in 1-D, Bessel function Jo(k\x—x'\) in 2-D, and 
sinc(fc|x— x'\) in 3-D. 

In Ref. |18| the interaction with electromagnetic fluctuations in metal has been considered. It has been assumed 
that each mode is characterized by an ohmic response. In our notations it corresponds to J(k,uj) = rj(k)-uj. Due to 
this assumption, the noise auto-correlation function becomes factorized at high temperatures, as in the DLD model. 
In particular, for the noise functional ( 5.1 0| ) one obtains 



•>N\r\ 



o PJo 



7T 



V (k) (sm 2 (kn-r(T))) Q 



(A.5) 



For the discussion of dephasing, we write again this result using the notations of Ref. |18| . It is assumed that 
77(A)) = e 2 /(trfc 2 ), where e is the charge of the electron and a is the conductivity. Integration over k leads to the result 



exp 



e 2 k B T 
2h 2 a 



\x a (r) ~ x b {r)\ 



2-d 



dr 



for d < 2 



(A.6) 



This result, that has been obtained in Ref. jl8| (equation (7.8) there), is very similar to the corresponding result 
(6.13) for the BM-model, the difference being the power (2— d) which is not universal. 

An apparently simpler example for an interaction with external bath that consists of extended field modes, is 
the electron-phonon interaction. Here we want to question the applicability of the BM model as an approximated 
description. The electron is assumed to be confined to a one dimensional "quantum wire" , while the phonons dwell 
in the 3-D bulk. Considering longitudinal modes, the coupling of the a — > (^1,^2,^3) "oscillator" with the electron 



is ki ■ x, with a coupling constant C a oc \Jk\ + k\ + fc| . Summing over a, as defined in ( 2.11 ), one obtains effectively 
J(uj) oc u) d+2 where d=3 is the dimensionality of the space where the phonons dwell. Thus a phonons-bath is similar to 
a superohmic BM-modcl. However, this derivation is somewhat miss-leading, since a cutoff uj c should be introduced, 
while for the original electron-phonon interaction a natural cutoff \x—x'\/c exists. 



B. AN INTERFERENCE GEDANKEN EXPERIMENT 



In this appendix we consider the interference phenomenon from two different points of view, and demonstrate their 
consistency. First we consider the free propa gati on of two-wavepackets superposition. The decay of the interference 



pattern will be dictated by the propagator (7.8). Then we consider the scattering of a simple wavepacket from a 
double-barrier. The suppression of interference paths will be dictated by the noise functional (Sec.VI(C)). Finally, we 
argue that both points of view are physically equivalent, and must lead to the same result, which is indeed the case. 

Consider a superposition of two Gaussian wavepackets which have the same momentum Pq, the same initial spatial 
spreading <7q, while their initial locations satisfies |i?02— R01 \— d. The Wigner function for this preparation is easily 
computed, and is of the form 

Pt ^{R,P) = ^G(R-R Q1 ,P-P ) + Ig(R-R 01> P-P ) + cos(^j^-Yg(R~(R 01 +R 02 ),P-P ) , 
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where G{x,p) denotes "minimum-uncertainty" Gaussian distribution, and 6P c =h/d. For free propagation p t (R, P) 
will develop fringes on the spatial scale Ax—^^-t . Note that this gedanken experiment is formally equivalent to the 
usual "two slit" diffraction experiment upon the definition AO = ^| = -p^/d. For propagation in noisy non-disordered 
environment the interference pattern is smeared on scale SP ~ v-t, due to the diffusive momentum spreading. The 
smearing factor is exp[— ^(5P(t)/SP c ) 2 ] leading to an exponential decay exp[— jp-t], that depends on the separation d. 

On the other hand, for propagation in noisy disordered environment, using (7.8), the exponential decay is exp[— ^rt], 
independent of geometry. 

Consider now scattering problem in one dimension. The potential is assumed to be V(x)=oo for x < 0, and 
V(x) = 6(x— 5) for < x. A simple Gaussian wavepacket with momentum —Pq and spatial spreading 00 <C d is 
launched from R . The scattered particle is detected at the range R which is assumed to be much larger than R . 
In the absence of noise and dissipation, it is not difficult to work out the explicit solution of this scattering problem. 
A "train" of Gaussian-like wavepackets will emerge from the scattering region. The spatial separation between each 
two wavepackets is d, and the probability density function will contain interference pattern in between. From the FV 
path-integral point of view (Sec.VI(C)), the interference pattern is due to the existence of "interference paths". In the 
presence of noisy disordered environment the contribution of these interference paths to the propagator is suppressed 
exponentially in time. 

Assuming that the dephasing during the time interval t sca ttering=Ro/v is negligible, one may use the results of 
the first gedanken experiment, leading to the same exponential decay. Thus, the result ( |7.8| ) for the propagator 
is consistent with the analysis of interference in Sec.VI(C). Our gedanken experiment can be used also in order to 
illustrate and clarify the observation that friction does not affect the interference phenomenon. This observation holds 
for disordered dynamical environment. 
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